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ABSTRACT 



X 

We propose a discrete model whose continuum limit reproduces the string sus- 
ceptibility and the scaling dimensions of (2, 4m)-minimal superconformal models 
coupled to 2Z)-supergravity. The basic assumption in our presentation is a set 
of super-Virasoro constraints imposed on the partition function. We recover the 
Neveu-Schwarz and Ramond sectors of the theory, and we are also able to evalu- 
ate all planar loop correlation functions in the continuum limit. We find evidence 
to identify the integrable hierarchy of non-linear equations describing the double 
scaling limit as a supersymmetric generalization of KP studied by Rabin. 
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1. INTRODUCTION 



The discrete formulation of N — 1 Superconformal Field Theories [1] coupled 
to world-sheet supergravity lags far behind its purely bosonic counterpart. Some 
of the results in [2] were obtained previously in terms of discrete models of 2D- 
gravity (see for example [3]). The supersymmetric extension of the analysis in 
[2] carried out in [4], [5] has no analogues in terms of random (super)-surfaces. 
In the double scaling limit [6], the theory of the KP (Kadomtsev-Petviashvili) 
hierarchy was shown to play a central role [7]. Motivated by this connection, the 
authors in [8] proposed an approach to the double scaling limit of 2D-supergravity 
coupled to superconformal matter using a supersymmmetric generalization of the 
KP hierarchy due to Manin and Radul [9] . An important feature of the one-matrix 
model wich will be central in our arguments is the fact that its partition function 
satisfies a set of Virasoro constraints [10], [11]. 

In this paper we propose a discrete model of 2D-supergravity with supercon- 
formal matter following the more "phenomenological" approach in [12]. Reasoning 
by analogy with the Virasoro constraints we define a discrete analogue of the 
Hermitean one-matrix model. Our basic postulate is to begin with a set of super- 
Virasoro constraints in the Neveu-Schwarz (NS) sector satisfied by the partition 
function. From them we derive the explicit form of our model and a set of super- 
loop equations. In the planar limit these equations can be solved exactly. In this 
way we compute the spectrum of anomalous dimensions which coincide with the 
(super)-gravitationally dressed dimensions of the (2,4m) minimal superconformal 
models in the NS and the Ramond (R) sectors. We compute arbitrary multiloop 
correlation functions on spherical topologies, and our results agree with those ob- 
tained in the continuum limit [13] using the super-Liouville formulation of the 
problem. Since there are no higher point functions computed in the continuum 
we cannot compare our results further. In spite of these encouraging properties, 
a derivation of our model in terms of "triangulated super-surfaces", orthogonal 
polynomials and generalized matrices is still lacking. Thus, the identification of 
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our model as a discrete version of 2D-supergravity should be taken as preliminary. 

We take as guiding principle in our work a set of super- Virasoro constraints 
satisfied by the partition function. One reason why we believe this to be a correct 
starting point is the prominent role the Virasoro constraints play in the description 
of the geometry of the moduli space of stable Riemann surfaces after Witten's 
work [14] on 2D-gravity and matrix models, and the proof of Witten's conjecture 
by Kontsevich [15] (see also [16]). According to Witten's theory, the intersection 
theory of certain line bundles on the moduli space M. g , n of genus g surfaces with 
n distinguished points is captured by the Virasoro constraints satisfied by the 
partition function. Any discrete version of N — 1 supergravity on the world-sheet 
should necessarily have to address similar issues for super-surfaces. In this case, 
however, the mathematical theory is not sufficiently well developed to allow us to 
borrow results which could shed light on the problem. We take a more radical point 
of view in expecting many properties of super-moduli spaces to be captured by the 
super- Virasoro constraints G n _iZ = 0, n = 0,l,... . S ince \G<m G m ^ ~ L n ^- m , the 
Virasoro constraints are automatically satisfied. The results of this paper can be 
interpreted as giving support to the validity of this basic assumption. 

We also find evidence that the generalized KP hierarchy appearing in our model 
is not the super-KP hierarchy of Manin-Radul [9] but rather a super-hierarchy 
defined by Rabin [17]. The basic difference between the two is related to the fact 
that the latter does not admit a simple presentation in terms of a Lax pair. We 
will comment on these issues at the end of the paper. 

This paper is organized as follows. In section two we collect several results 
concerning matrix models, loop equations and Virasoro constraints and present 
them in a way which will be easily generalized later. In section three we study 
the super- Virasoro constraints, we use them to derive the exact form of our model 
and then we derive the super-loop equation in its planar approximation. Section 
four analyzes the solution to the planar loop equations, the scaling limit and the 
spectrum of scaling operators. We obtain the dressed gravitational dimensions 
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expected in the Neveu-Schwarz and Ramond sectors of the (2,4m) minimal super- 
conformal models coupled to 2D-supergravity. As an appplication we compute in 
section five correlators of arbitrary numbers of loops (bosonic and fermionic) in 
planar topologies. Section six contains our conclusions and out-look, and our re- 
marks on the connection between our model and the super-KP hierarchy described 
by Rabin [17]. 

2. VIRASORO CONSTRAINTS AND LOOP EQUATIONS 

We review succintly in this section some properties of the one-matrix models 
and their loop equations. The starting point of Kazakov's analysis of multicritical 
points [12] was the planar loop equation* 



E fc Jfcajfc=T w (0 = / dl ' w{l - z >(0 C 2 - 1 ) 
k>\ 



w(l) describes a loop of length / bounding a surface with the topology of a disk. 
This equation follows from some simple heuristics, but it can also be derived from 
a Hermitean matrix model [18]. We take the partition function to be 

/A T 
d Na $exp[-—trV(Q)] 

V($)=Y,9k* k A = e^ B (2.2) 

fc>0 

where $ is a Hermitean N x N matrix and hb is the bare cosmo logical constant. 
The loop operator is represented by 

A in a in 

w (l) = -tre 1 * = Y --tr§ n = V -w^ (2.3) 

n=0 n=0 

Writing the partition function in terms of the free energy Z = e N F , F = Fq + 



* In the sequel we will use the same symbol w(l) for the loop operator and its expectation 
value. From the context it will be clear which is the correct interpretation. 
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N 2 Fi + N 4 i<2 + the (expectation values) moments are given by 

«,(«)= A W ^ = -A 2 ^- (2.4) 

dg n 

Near the critical point /i c , the genus expansion of (2.2) behaves according to 

Z0m b ) = Zh[N 2 (^B ~ /ic) 2 " 7 "] 1 "* (2.5) 
h 

where h is the handle counting parameter and 7^ is the string susceptibility. Intro- 
ducing the renormalized cosmological constant /is — /i c = a 2 t, where a is a cut-off 
with units of length, the double scaling limit [6] is obtained by taking N — > oo, 
a — > 0, and keeping fixed the combination 

Na 2 ~ rst = k _1 (2.6) 

which is the string coupling constant. (More details and references to the literature 
can be found in the review articles [19], [20]). For later convenience we derive the 
planar-loop equations (2.1) through the Virasoro constraints satisfied by (2.2) [21]. 
They are obtained by making the change of variables $ — > $ + e4> n+1 . After some 
simple manipulations we obtain 

L n Z = n > -1 (2.7) 

with 

A 2 n d 2 d 

L " = ^g^^ + g* 9i ^ (2 - 8) 

To leading order in (the planar limit) we write Z = e N2p{) and obtain 

2 -A<9F 8Fq 8Fq 

A 1^1T~ a + l^ k 9kj; = 2.9 

~ ogk og n - k ^ dg k+n 

Using (2.4) it is easy to see that (2.9) is exactly identical to (2.1) . Hence the planar 
limit of the Virasoro constraint (2.7) is nothing but the planar loop equation. The 
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equations (2.1) (2.9) are solved by introducing the Laplace transform of the loop 
operator 

00 



w 



e-^(0* = E (2 - 10) 




with the assumption that behaves well at / = and 00. With the definition 
(2.10), (2.1) becomes an algebraic equation 



w{p) 2 - V'(p)w(p) + Q{p) = 



fc-i 



fc>l r=l 

Two remarks should be made at this point. First, if we make the identification 

AV2 d N . 

a n = tttt- , n > ; a~ n = -—^=ng n , n > (2.12) 

N dg n AV2 

the L n 's can be rewritten as the components of the energy- momentum tensor of a 
free massless scalar field 

Ln = lj2 :Ct -k a k+n- , T{p) = Y,Lnp- n ~ 2 (2.13) 



neZ 



Defining 



dX(p) = E ctnV^ 1 = dX+ + dX- 



neZ 



dX+ = J2*n P - n - 1 



(2.14) 



n>0 



we obtain 



(n) 



lim ^-Z^dX+Z oc w(p) = V 
lim -J-^dX^Z oc V'(p) = V kg k p k ~ x 



(2.15) 
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and the Laplace transformed equation (2.11) becomes 

in the limit p — > 0. This is important because with the identification (2.15) we can 
find the potential for the matrix model. The planar loop equations are therefore 
equivalent to (2.16) . The second remark has to do with the form of Z in terms of 
eigenvalues. If the eigenvalues of $ are Ai, Ajv, (2.2) can be written as [22] 

/AT 
n ^ - x ^ 2 ex p[- a e v ^ ( 2 - i? ) 



We could instead write 



/AT 
Hd\iA 2 (\)exp[-— ^^(Ai)] (2.18) 

i 

and use the Virasoro constraints (2.7) (2.8) to determine the form of A. Indeed, 
if we act with (2.8) in (2.18) and perform some simple integrations by parts, we 
obtain a set of differential equations satisfied by A 

whose solution up to an irrelevant constant is A = Hi<j(^« — x j)- The solution to 
the planar equations (2.11) takes the form [12] 

w(p) = \ (V(P) - M(p 2 Wp*-R) (2.20) 
in the case when the potential is even V{p) = V(—p). Since 

w(p) = ^ + ^ + ... , «,<°> = A(i2) (2.21) 
p p 6 

M(p 2 ) is completely determined by requiring the right hand side of (2.20) to have 
only negative powers of p in the large p expansion, and R is completely determined 
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in terms of A. Kazakov showed [12] that 

dw {p) 1 



(2.22) 



<9A Vp^R 

A = e-»° = A(R) = £ kg 2k R k ( ^ ) ^ = / dxV (* + ^) ( 2 - 23 ) 

At the m-th critical point 

A = 1 - (1 - R) m = e-» B = l-a 2 t 

(2.24) 

1-R = a 2 > m u 

and t and u are scaling variables. The variable t is the renormalized cosmological 
constant and u is the "heat capacity" of the theory. 

To compute loop correlators notice from (2.22) that 

Taking the limit as k — > oo, and defining the renormalized length according to 

I = ka 2 ' m (2.26) 

(I kept fixed as a — > 0, k — > oo), using the scaling limit (2.6) with 7 s t = — ^ and 
using the scaling variables (2.24) we obtain 



oo 



iu(Z) = y/n(tr$ 2k ) = / dt'e~ lu ^ (2.27) 



To compute multiloop correlation functions all we need to know is 

MZi)w(Z 2 )) = 7r(tr$ 2fcl tr$ 2fc2 ) 

A /- 9 I II \\ 7 7 -2/m ( 2 - 28 ) 
A 002fc, 

Perturbing the string equation (2.24) with the g 2 k 2 coupling and taking k 2 — > oo 
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as in (2.26) we obtain 



dg2k V vr 

which together with (2.27) , (2.28) yields 



- = -o- 2 " 1 /" A-'»« (2.29) 



(w(h)w(l 2 )) = Vhh e -r— t- (2.30) 

n + h 

Repeating the same arguments we obtain 

(W(h) . . . W(l n )w(l n+1 )) = -Ky/l^e- l ^ U ^- t (w(h) . . . W(l n )) (2.31) 

in agreement with the results in [23]. 

In the generalization to the supersymmetric case in later sections we will follow 
closely the arguments in this section. 

3. SUPER- VIRASORO CONSTRAINTS 
AND SUPERLOOP EQUATIONS 

We now introduce the analogue of w(l) which we take to depend on two vari- 
ables, w(l,9), a bosonic and a fermionic length I and 9 respectively. We can imagine 
these two paramenters as characterizing the boundary of a super-disk. As in the 
previous case we can introduce the super-Laplace transform 

oo 

w(p, II) = v(p) + Uu(p) = J dl J d9e- pl - nd w(l, 9) (3.1) 



Some properties of the Laplace transform which are useful are the following: 

oo 

C[f] = f(p, H) = JdzJ d9e-r z -™ (f (z) + 9h(z)) 
o 

= /i(p) + n/ (p) 

oo 

fi(p) = J dze^Mz) , i = 0,l (3.2a) 



-u(h+l 2 ) 
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£p/] = (n+„!j)£[/]-/o(0) , D = ^ +0 l. (3.26) 



:|r+n|)£[/] (3.20 



£[/»<,] = -(-) a /£[/]£[ g ] (3.2d) 

z 

(/ o ^ e) = y ^ y ^7(^', e') g(z -z',9- e') 



= h ° 50 - /o o gi - o gi ) (3.2e) 

where (9/ is the grading of / (df — if / is even and (9/ = 1 if it is odd). In 
the last line fi o ^ is the standard convolution of functions. We introduce also the 
symbol 

P = «-*| (3-3) 

related to £) = ^ + via the Laplace transform. 

The first (of three) derivation of the planar loop equations is based on the 
analogy with the c = 1 energy-momentum tensor. Assuming again that the loop 
w(l,6) behaves well as I is near or oo, we can expand the Laplace transform 
w(p, IT) in inverse powers of p 

v — -\ v — -\ 

fe>0 V k>0 y 

Consider a c = 1 free massless superfield X(p, U) = X(p) + Uip(p). Its super- 
energy-momentum tensor is 

Tip, n) oc: DXdX := ipd p X + IT : (d p Xd p X + 6^) : (3.5) 

Writing 

0X(p) = 5>«p- n - 1 , ^(p)= £ 6,p- r - 1/2 (3.6) 

neZ reZ+1/2 
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we can identify a n , b r with bosonic and fermionic couplings 

A d N n , n , n „, 

ap = ~N~dg~ p ' a ~ p = ~li P9p ' P = ' ' ' ' ' ' ( ^ 

A d N 
bp+1/2 = ' & -p-i/2 = > P = 0,1,2,... (3.8) 

The Laplace transformed loop w(p, U) is identified with the positive frequency part 
of w(p, n) ~ DX + and the analogue of the potential in the one-matrix model is 
identified with DX~, more precisely DV(p : II) ~ DX~. This leads to 

v( P , n) = + £ fc+1/2 n/) (3.9) 

fc>0 

Writing Z = e N F we can also identify the moments vS- k \ with derivatives of 
F: 

u< > = A , = -A 2 |^ ; = -A 2 -^- (3.10) 

dgn din+l/2 

With these identifications, we can take the planar limit as in the pure gravity case 

lim N- 2 Z- 1 T(p,U)Z = (3.11) 

with Z ~ e N F °. Some simple algebra yields what should be considered as the 
Laplace transform of the superloop equations: 

KDw + DKw = wDw + Q 

K = DV = J2(Zk-i/2 + Kkg k )p k ~ l , Q = Q 1 + UQ 
k>\ 

fc-1 

^ = 2^^{W" i_1) + (k- i-jm k _ 1/2 v^-vy- 1 

(3.12) 

= EE^ N " 1) +^i / /- H) }f 1 
fc>i j=i 

Using the Laplace transform formulae in (3.2) one shows that (3.12) is the Laplace 
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transform of the superloop equations 



VKw(l, 9) + 21CVw(l, 6) = (wo Vw)(l, 6) 

K^^gkde + ^k+i^df- 1 (3.13) 
fc>l 

which is strongly reminiscent of (2.1) the starting point of Kazakov's analysis. 
It should be quite interesting to derive (3.13) heuristically in terms of gluing su- 
perdisks through their boundaries. 

In the second and more fundamental derivation of the superloop equations 
(3.12), (3.13) we begin by constructing an "eigenvalue model" similar to (2.18) 
and use super- Virasoro to determine the measure. The form of the potential (3.9) 
suggests the introduction of N pairs of eigenvalues (Aj, 0$), one bosonic and the 
other fermionic. The potential for this eigenvalue model is taken to be 

V(X, 9) = J2 + &+i/ 2 0iAf) (3.14) 

fc>0 i 

and the partition function is written as 

/AT 
II d\id9iA(\, 6) exp[--y(A, &)} (3.15) 

i 

The explicit form of the super- Virasoro generators in the 6=1 case with the 
oscillators (3.7) , (3.8) is given by 

Or = ^ E ^ + E^ + g kn ^- (3,6) 



A 2 n d 2 °° d 

L - = VN2 E x m x n _ , + E k 9k 



2N2 ^ dg k dg n - k ~ dg k+n 



A 2 n ^i 2 ' n d d ( n d 

r—1/2 r=l/2 



(3.17) 



Since L n is obtained in terms of anticommutators of GV's, it suffices to impose on 
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Z only the fermionic constraints G iZ = 0. It is convenient to write G 1 as 



t^ d 9k+n ^ 0&+n-l/2 ^ ^ 9^+1/2 fyn-l-fc 

If we recall the explicit representation of the action of the algebra of supervector 
fields (super- Virasoro without central extension) on the space of functions of (A, 9): 

9„-l/2 = -0\ n §^ + A"^ (3.19) 



we find that the action of G n _i/2 on Z can be traded off by the action of 

i 

on the exponential term. Integrating by parts and using identities like 

n—l yn yn 

we obtain a set of differential constraints on A 

whose unique solution (up to an irrelevant multiplicative constant) is 

A = - Xj - Ofy) (3.24) 

i<j 
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Hence our model is explicitely given by 

/AT 
Yl dXidOi J](A. t - \j - 9i9j) exp[--F(A, 9)\ (3.25) 

i i<j 

Introduce now the explicit representation of the superloop operator 

Hi,e)^^J2 elK+m ( 3 - 26 ) 

i 

and its expectation value with respect to Z. Acting on (w(l,9)) with the oper- 
ator (VIC + 2JCV) appearing in (3.13) , using the super- Virasoro constraints and 
the factorization of amplitudes in the large N limit we obtain (3.13) after some 
computation. 

Finally we could obtain (3.12) , (3.13) by using (3.10) , and the explicit formula 
for Z~ x G n _xZ = 0, n > 0, and Z~ x L n Z = 0, n > -1, with Z ~ e N2p ° and taking 
the large iV limit. This gives the third derivation of the planar loop equations 
(3.12) which we can write in components according to 



a) («(p) - V\p)) 2 + (v(p) - H(p))'(v(p) - ^(p)) = A 

b) (v(p)-ap))(u( P )-V'(p)) = A 1 



(3.27) 



where 



k>0 

V\p) = Y J k9kP k - 1 



(3.28) 



k>l 



and 



(3.29) 



A = V'(p) 2 + Z\p)ap)-Qo 
A 1 = ap)V'(p)~Qi 

Qq and Qi are given in (3.12) although their explicit form is unnecessary except 
for the fact that they are polynomials in p. The primes denote differentiation with 
respect to p. 
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After the three derivations presented of (3.27) the next step is to solve the 
superloop equations, take the scaling limit and find the critical points and critical 
dimensions. This we do in the next section. 

4. SOLVING THE SUPERLOOP 
EQUATIONS. SPECTRUM OF THE MODEL 

In solving the superloop equations we will make the simplifying assumption 
that the bosonic part of the potential is even 

0) = Y, E^fcAf + &+i/2*iAf) (4.1) 

fc>0 i 

Before attempting the solution of (3.27) there are a number of useful remarks 
that should be made concerning the solution of the 2_D-gravity case (2.11), (2.20). 
As written in (2.11) it seems that the polynomial Q contains a number of "initial 
conditions" for some of the loop moments w^ k \ If the potential is of order k, 
apparently the first k — 1 moments are completely arbitrary, and all other moments 
are computed in terms of them. 

This is an incorrect impression if we think in terms of the original matrix 
model. Once the couplings and the cosmological constant A are given there 
are no ambiguities. If we think in terms of the formal perturbative evaluation, we 
expand the potential in powers of all couplings gk with k ^ 2, i.e. we keep the 
quadratic term g2& 2 in the exponent and evaluate the free energy F using Wick 
contractions. The different correlators are then formally analytic in all couplings 
gk , k 2, and there are no free constants. This can also be shown by solving the 
Virasoro constraints perturbatively in g^ , k ^ 2. This uniqueness of the formal 
loop correlators when written in terms of (gk, A) implies that the solution to (2.11) 
is unique provided: i) it is parametrized by A and it is perturbative in the couplings 
gk , k 7^ 2; ii) w(p) only contains negative powers of p for \p\ large; and of course 
iii) w 2 (p) — V'(p)w(p) is a polynomial, which is another way of expressing (2.11) . 
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Thus although different ways of parametrizing the cut in (2.20) and the function 
M(p 2 ) might seem to yield inequivalent answers, when we express R in terms of A 
through the condition 



w(p) = - + 0(\) , A = A(R) (4.2) 
p p z 



all the answers will be the same. 

We can take the previous remarks and apply them in our situation. It is again 
true that the perturbative evaluation of Z, or the perturbative solution to the 
super-Virasoro constraints yields unambiguous answers for the superloop moments 
tS- k \ v^ k \ Thus we stress the fact that if we find a superloop w(p, II) = v(p)+Uu(p) 
satisfying: i) at large \p\ it only contains inverse powers of p; ii) it is perturbative 
in all couplings except the quadratic even coupling; iii) the left hand sides of (3.27) 
computed with the proposed solution w(p, II) should be polynomials (Ao and Ai); 
then when we use the string equation 



«(p) = - + 0(l) (4.3) 
p p z 



to express the auxiliary parameters in terms of A, the solution is unique. This 
uniqueness property of w(p, U) is very important and it will be used presently. 
The solution may be parametrized differently but the answers in terms of A, 
k 7^ 2, £fc will be the same. 

To obtain a preliminay idea of the possible analytic structure of u(p), v(p), we 
solve (3.27) explicitely in terms of Aq, Ai. The second equation in (3.27) gives 



" w - {W = SofHo (4 ' 4) 

which turns the first equation in (3.27) into a quartic equation in u — v' . However 
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since Ai is odd, A^ = 0, we obtain after some simple manipulations 

nip) - V (p) = - ^ (4.5) 

«(p) - ^(P) = 4f= (4-6) 
v A o 

In principle Ao, Ai are respectively even and odd polynomials in the fermionic 
couplings and a similar conclusion applies to u{p) and v(p). Among the 

two non-trivial solutions to the quartic equation defining u — V we choose the 
one which gives a contribution to u — V' of order zero in the fermionic variables. 
Imitating the solution to the planar bosonic case we choose u(p) to have a single 
cut 

u(p) = V'(p) - M(p 2 ) v / p 2 -R + ... (4.7) 
and for the fermionic operator v{p) we take the starting ansatz 

v(p) = UP 2 )+UP) - N^p 2 )y^R- + • • • 

VP ~ R (4.8) 

e + (p 2 ) = E^+l/2P 2fc , UP) = J2^ + 3/2P 2k+1 

k>0 k>0 

From the form of (4.5) we easily guess that u{p) must also have a contribution 
proportional to (p 2 — R)~^ . The polynomials M, N± are chosen to cancel the 
positive powers in g, £+, Substituting (4.8) into the left hand side of (3.27b) 
we see that we find a polynomial on the right hand side. As expected, with (4.7) 
we find that we need to modify the ansatz. It is easy to check that it suffices to 
add to u{p) the piece 

R n+(R)(n (R)+pn+m (4Q) 



M(R) (p2 _ R )3/2 

to guarantee that the left hand side of (3.27a) is a polynomial. In (4.9) n^(R), 
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n 



i(R) and M(R) are obtained by expanding N±, M in powers of p 2 about p 2 = R: 



M(R) = M(p 2 )\ p2=R 
N ± (p 2 ) = n±(R) + nf(R)(p 2 - R) + . . . 

To summarize, our solution to the planar superloop equations is given by 



(4.10) 



W = V(p) - m^Vp^R- ^ 



v(p) = ZAp 2 ) + Up) - nAp 2 )Vp^R- ^== (4-12) 

VP - R 

As we will see below, N± are linear in the fermionic couplings and u(p) is at most 
bilinear in fermions. This is rather surprising and one might be tempted to add 
higher order terms in the ^-couplings. Note however that (as soon as we give 
the explicit form of M, N±) u(p) , v{p) have only negative powers of p as \p\ — > 
oo, that the solution is perturbative in the couplings k ^ 2 and furthermore, 
by construction the left hand side of (3.27) are polynomials. If the reader has 
accepted the uniqueness arguments given at the beginning of this section he or she 
is unavoidably led to conclude that after R is traded by the cosmolgical constant the 
same conclusion holds. Therefore any other (more complicated) dependence in the 
fermionic couplings is spurious and it will be redefined away when the correlators 
are expressed in terms of A. This is purely a planar phenomenon and it is expected 
not to hold to higher orders in the large N expansion in our model (3.25) . We 
will offer a more geometrical explanation of this fact in the next section, after we 
evaluate correlation functions of loop operators. Nevertheless this feature of the 
planar superloop equations is rather surprising. 

The determination of M(p 2 ) is exactly as in the bosonic case, and this implies 
therefore that the critical points will be labelled by a positive integer m with a 
string susceptibility 7 s t = — Decompose u(p) into uo(p) + U2(p), uo(p) (resp. 
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U2(p)) is of order zero (resp. two) in the fermionic couplings 



MP) = y'ip) - M(p 2 )^p>-R (4.13) 



This form of uq(j)) automatically implies the expansion 



, , A(R) 

u (p) = ^ + _ + . .. + — + ... (4.14) 



We can arrive at the same result directly from our model (when V(p) is even) 
by showing that all odd expectation values (J2i^i k+1 ) vanish when all fermionic 
couplings are set to zero. In this case, the partition function of our model takes 
the form 

r 2N n 

Z= I Y[d\ l d9 l l[(\ l -\ j -9 l 9 J )exp[--Y,Yl92k\l k } (4.15) 

i=l i<j k i 

The integration over #'s yields A(A)P/(A- 1 ), where A is the standard Vander- 
monde determinant and Pf(X~j 1 ) is the Pfafhan of the antisymmetric matrix 
Mij = (Xi - Aj) _1 i ^ j, Ma = 0. Since A(A)P/(A ij - 1 ) is even under A, -> -A, it 
is clear that any correlators of the form X 2k+1 ) vanishes. 

Multiplying (4.13) by p and introducing the new variable t = p~ 2 we have 

F(t) = t-^uoit- 1 ' 2 ) = g^ 1 ) - M(r 1 )r 1 v / T^itt 
g(p) = pV'(p) 

The left hand side F(t) is analytic near t = 0: F(t) = A(R) + u^ 2 h + . . .. Similarly, 
for t small F(t)(l — Rt)~% is also analytic near t — 0, hence — Rt)~* — 
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M{t 1 )t 1 should be analytic as well and we can determine M 

M(t- l ) = Y,M n r n , 



9^) lA ^ 



We need in particular 



fc=0 

Furthermore 



^ , 2k\ 2k + 2 h 

M ° = E J -$-&9m-2 - /(*) 



by analyticity the term in parenthesis in (4.18) is independent of t and can be 
computed using f(R). From (4.14) , (4.16) we conclude 



dF(t,R) _ d R A(R) 
dR ' " ^/l^Rl 



(4.19) 



9rMR) = \f{R) + (4-20) 

H k f 2k\ r r 
MR) = E k 92kjr , H + 4-) (4-21) 

fc>i V K / o 

Evaluating the term in parenthesis in (4.18) (which is independent of t) at t^ 1 = R, 
we obtain the useful identity 

M(p% 2=R = 2d R A(R) (4.22) 
The determination of the fermionic functions N±(p 2 ) goes along the same lines. 
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Defining 



„ ^ v (2fc+l) v (2k) 



fc>0 y fc>0 y 



(4.12) splits into two equations 



v+(p 2 )=Z+- , v-(p) = e- - N-y/^R (4.24) 

vr - R 



The equation for V-{p) is identical to the bosonic case. Introducing t — p 



-2 



v+(t- A ) = - (4.25) 

since Vl — Rtv + (t^ 1 ) is analytic near t = 0, the analogue of (4.16) becomes 

y - Rtt l - l dt = j> N+it-^t^dt (4.26) 



and 

As v+(t, R) is divisible by t, the term in parenthesis is independent of t, and only 
the zeroth order in N + = J2 n =o N+ n t~ n is relevant 

MR) - -\n +0 + ^ = -iyj w f (") («8) 

fc>o V h J 

Expanding N + about p 2 = R 

N+(p 2 ) = n+(R) + n+(R)(p 2 - R) + . . . (4.29) 
and using the t-independence of the term in parenthesis in (4.27) we obtain 

n+(R) = -2p' + (R) , n+(R) = -4p'l(R) (4.30) 
We can carry out similar computations for V-(p), and determine as well U2{p). The 
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final results can be summarized as follows: 
du (p) 



dA ^ 



i r>k / 2k\ 

W=R ' m %^(k) (431) 



dv+(p) PP'+(R) n>(m _ R h (Zk\ 



Expanding in powers of |, we can write equivalently 

M 2 2fe) = 4l>i*>2 = K ) p'+W-(R)dAR , k > 1 (4.38) 



= = - 4 ^ + (j PWWRdAR (4.39) 

The subscripts or 2 in (. . .) indicate orders or 2 respectively with respect to the 
fermionic couplings. Equations (4.35) - (4.39) are the basis for the computation of 
loop correlators in section five. 
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To take the scaling limit of the previous expressions we notice that the zeroth 
order equations in fermionic couplings are those of the purely bosonic theory. Fur- 
thermore A(R) in (4.31) is independent of any fermionic couplings. In analogy 
with the analysis in [12], near the m-th critical point A(R) is given by 

A(R) = 1-(1-Rr + J2tn^-R) n 

(4.40) 

= A = e^ s ~ 1 - a 2 t 

where t is the renormalized cosmological constant and is the bare coupling of the 
n-th scaling operator. Exactly at the m-th critical point all t% — 0. Introducing 
the scaling variable 

R = l- a 2 ' m u (4.41) 
The string equation and the renormalized couplings are as in the bosonic case 

u m -^2t n u n = t (4.42) 



tl = a 2 ^- n / m h n (4.43) 
with string susceptibility at the m-th critical point given by 

7 st = -1/m (4.44) 

and in the double scaling limit we again keep 

Na 2+i/m = K -i ( 4>45 ) 

fixed. From (4.35) and (4.40) we can determine the form of the bare scaling bosonic 
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operators following Kazakov [12]: 



2(n + l)A^ ( 2k \ v , i=1 



There are two contributions to the planar expectation value: 

(0 = <Oo+<0 2 (4.47) 

Notice that the contribution to A(R) of adding the operator (4.46) to the potential 
is 

5A(R) = — d R (l - R) n+1 = (1 - R) n - (1 - R) n+1 (4.48) 

n + 1 

and in the continuum limit only the first term survives. This agrees with the 
identification of scaling operators in [24]. From (4.43) the renormalized operator 
is 

<r+ = a 2{l - n ^a^ + (4.49) 
Since d/dA = —a~ 2 d/dt we obtain 

d '^° = -i^TT (450) 

and 

(<7+) ~ i 1 + 1 M+n/m = t 2- lst+ {d n -l) ^ ^ 

We can identify the gravitational dimension of <r+ as 

rfn = - (4.52) 
m 

This result coincides with the scaling dimensions in the NS-sector of a (2,4m) 
N = 1 superconformal minimal model coupled to 2_D-supergravity. 
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For the fermionic operators, in analogy with (4.40) we introduce the fermionic 
scaling variables according to 



dp± 



tI^^-RT (4-53) 



dR 1 2 

n 

This definition determines the form of the fermionic operators to be 

n ( A\k / n \ N 




'» B - = -tEtS44^-E^ (4.55) 



A£j(» + l)/», ,. = | 



Hence 



^ (^(^n ± )) = p / T (^)5 A i?(l - RT (4.56) 

To unambiguously determine the scaling behavior of the fermionic variables 
we find it convenient to work first with the odd bosonic operators i.e. those de- 
termined by the odd expectation values ^A? fe+1 . From (4.39) we see that we 
have three point functions involving one odd bosonic operator and two fermionic 
operators u + . The explicit form of the odd bosonic operators is 




B _ (-4) fc \k, „ 2t+1 



To derive this expression we must go back to the original bosonic loop equation 
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(3.27a) and add a small odd perturbation to the bosonic potential: 

V'ip) = k 9kP k ^ = Vi( P ) + VL (p 2 ) (4.58) 



k>i 



VUp 2 ) = J2( 2k + 1 ^2k + iP 2k (4.59) 

k>0 

The loop equation has a solution (to first order in V!_): 

«o(p) = V[{p) + VUP 2 ) - M{p 2 )^p^~R- VM Jfl. + 0{V a ) (4.60) 

VP - R 

Following arguments similar to those explained previously, and defining A_(i?) by 

u(p) = A(fl) + A4R) + L 

P p £ p 6 

we obtain: 



dKU ° = ~^=t + ' *<-(PX'=n = -^- (4-62) 



and 



k>0 

Therefore 



««A- = ~E^T^( 2 t V (463) 



^(^E^)^C)^ (-) 

We see that the operator (4.57) behaves according to 

^ K (^(aS-))=dRA-(R)d A R(l-R) n (4.65) 
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d R A- = -l'£ t S~( 1 - R ) n ( 4 - 66 ) 



2 

n 



We can now derive the scaling dimension of by computing two-point functions: 



" " oij! 



^ (^K B -<->) = ^(1 " i?)^ A i? (4.67) 



Hence 



dt(Vn~Vp~) = l _ t Vm+nlm+p/m-l a 2{n/m-l) a 2{p/m-l) 

ztyik 



and 



S- 2(1-^) 



(4.69) 



Therefore the dimensions of the odd operators can be read off from 

(<7-<7p) ~ iVm+n/m+p/m ( 4J0 ) 

coinciding with the dimensions of the even operators (4.52). This duplication of 
operators is familiar from the bosonic case. 

The bosonic loop u(p) in our solution contains a term bilinear in p+. Equation 
(4.39) together with (4.57) implies: 



N 2 

(*S-) = -^R P '+(R)p'i(R)dAR(l -R)' 



k 

(4.71) 
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From the scaling behavior of p'± (4.53) we obtain 



at2 

(^~^ + ^ + ) = jj(q - p)(i - RY +q+n - l d^R 



hence 



/ b— b+ b+\ _ (q-p) f(l/m—l)+n/m+p/m+q/m—l/rn2(n/m+p/rn+q/r 



Since for scaling operators we expect 

(JJOi)^* 2 -^^*- 1 ) 



we conclude 



and 



We therefore identify 



d+ + d+ = -(p + q-1) 
F y m 



p m 2m 



Finally, from (4.36,37) and (4.54,55): 

^(^-)) = -// + («l-fl)« 



implying 

QJjyB+^B-x = 1 f fl/m-l)+p/m+q/m a 2(p/m-l)+2fq/m-l) 

Thus 
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and 



Summarizing 



and these are the gravitational scaling dimensions of the operators in the Ramond 
sector for the (2, 4m)-minimal superconformal model coupled to 2_D-supergravity. 
Notice that in the spectrum (4.82) all states are doubled except for the state with 
dimension cLq. This is a state in the boundary of the Kac table analogous to 
the redundant operator a m -\ in the m-th critical bosonic model. It is tempting 
to interpret u Q as related to the ground state of the theory. We also obtain the 
precise scaling behavior of p' ± : 

P± = T^^r^u) (4.83) 
r ± («)=X;r n ± « n , r^r^a-^-^ (4.84) 

n 

The reason why we had to go through such a long argument to compute the 
fermionic scaling dimensions is related to the extra term ±l/2m in (4.82). If we 
had only considered (4.79) the identification of the Ramond sector would have been 
ambiguous. To distinguish between different possible assignments we had to study 
carefully the coupling of a~ to v+v+. We have therefore shown that our model 
in the planar limit has critical points labeled by m — 1, 2, 3, . . . and with scaling 
dimensions in agreement with the NS- and R-sectors of the (2,4m) superminimal 
models. It is also easy to see that in the NS-sector the one-, two- and three-point 
functions of scaling operators agree with the results obtained in the continuum 
super-Liouville theory [13]. We can compute arbitrary correlators as well. In the 
next section we compute the correlation functions for an arbitrary number of loop 
operators. 
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5. PLANAR LOOP AND SUPERLOOP 
CORRELATION FUNCTIONS 



As an application of the result in the previous section we compute correlation 
functions of loop operators. We start with bosonic loops. We rescale for conve- 
nience by a factor of Z -1//2 the definition of the loop operator in (2.27). This section 
is a translation of (2.25-31) to the present situation. Define 

r N 

We will take a — > and k — > oo keeping I fixed. There are two contributions in 
the planar case to (5.1) coming from (4.35) and (4.38) 

U(l) = U (l) + U 2 (l) (5.2) 

For Uq(1) the result is as given in (2.27) apart from a trivial rescaling 

oo 

(W)) = li jdt' e-W (5.3) 

t 

The other contribution comes from 

4 2fc) = ^ -^ P ' + {R)p'-{R)R k dKR (5.4) 

Using (4.41,83) we obtain 

(U 2 (l)) = -- T + (u)r-(u) e~ lu d t u (5.5) 

K 



In total 



oo 

(U(l)) = i- fdt' e~ lu ^ - -e~ lu dtu t+(u) t-(u) (5.6) 

Kl J K 
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Similarly we can introduce two fermion loop operators 

N I — N 

v+(i) = 2v^J>a?* vL(0 = \jJ2 9 ^ k+1 ( 5 - 7 ) 
t=l v t=l 

After some simple computations using (4.40,41,45,83) we arrive at 

oo 

(V±(Z)> = ±- / dt' e- lu{ t } r T (u) d t >u (5.8) 

« y 

It may seem surprising that (V±(Z)) contains a factor a 1 /" 1 . However, when we 
think of superloops with length I and (bare) superlengths 8±: 

W ± (l : 9) = U(l)+9iV ± (l) (5.9) 

we define renormalized superlengths 

Q ±= QB a l/m ( 51() ) 

and they behave dimensionally as Z 1 / 2 , as we might expect. It is possible to show 
that U (Z) and V±(l) admit an expansion in terms of microscopic scaling operators. 
The bosonic scaling operators (4.46) have expectation values given by (4.50) and 

^(vSh = -2[(1 - R) n - (i - RT +1 } P ' + (R)p'-(R)dAR 

as one can show using (4.38). In the continuum limit at the m-th critical point: 

1 u n+1 

2k 2 n + 1 ( 511) 



dt{(T n )o 



This together with (5.6) yields 

U{1) = 2k V(-l) n+1 — o n + (singular as Z -> 0) (5.12) 
The singular terms are analytic in t and they correspond to contributions from 
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microscopic loops. Similarly in the fermionic case: 



d A 



^(^=- P ' T (R)d A R(l-RT (5.13) 



In the continuum limit 



d t (v±)=±— 1 T^{u)d t u (5.14) 



'2k 2 

and comparing with (5.8) we obtain 



_ 00 (-) n + 1 

V±(l) = 2Ka 1/m , inv n + (singular terms as Z -> 0) (5.15) 

n=0 n ' 



Thus 



w±(z, = t/(z) + egv ± (i) = 2kJ2 i n {° n + e ± ut) (5.16) 

n=0 



This result is dimensionally consistent with the dimensions computed for a n ,u^, 
U ~ r 1 / 2 , V + ~ Z 1 / 2 , VL ~ r 1 / 2 . Then we use W± = U(l) + 6%V±(l) where 9 B is 
dimensionless but 9+ ~ Z _1 . Hence [0+] ~ Z~ 1//2 and ~ Z 1 / 2 . Since [Z] = — 1/m 
in units of t, [9±] = ±l/2m and 



dim z/^ = dim <r n =F ^ — (5-17) 



as expected. 

Finally we compute multiloop correlation functions. We need to use equation 
(2.29), and similar formulae for r ± . It is not difficult to show that when k — > oo 
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as in (2.26), 



A du 
Nd^k 




(5.18) 



A_?LL = Ka iMl^ (5 19) 



= 2Ka 1 /™./I e -'« (5.20) 

After some simple algebra one finds 

(tf(Zi) • • • U(l n )) = (-2^) n - 1 (C/(Zi + • • • + Z„)> (5.21) 



(V±(/i)f/(/ 2 ) • • • U(l n )) = {-2Kd t ) n - l ^{U{h + ■■■ + l n )) (5.22) 



(V+(h)V-(l 2 ) ■ ■ ■ U(Q) = (-2«)- 1 9r 3 ^^ r (C/(Z 1 + • • • + l n )) (5.23) 



From these expressions we can derive the correlator of n superloops with n + W + 
operators and rT W- operators. Define the differential operators 

Vt = d t + e ± ~ (5.24) 



Then 



(W+(h, 6{) ■ ■ ■ W + (l n+ 9n + )W-(l n++h e n++1 ) ■ ■ ■ W-(l n , e n )) = 



n+ n 

{-2K) n ~ l \{vj \{ VJ d -i (U( h + l 2 + ... + l n )) (5.25) 

i=l j=n+ + l 
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where 



00 




(5.26) 



t 



With this definition U(l) in (5.6) becomes 

(U(l)) = ^(-d^ + r + r-d t )e- lu (5.27) 

Notice also that 

n^-^+Mr 1 ^ (5-28) 

with 

n± 
i=l 

Therefore the superloop correlators (5.25) depend only on the total length L = 
l\ + I2 . . . + l n and the total superlengths 9± = &±i- A posteriori this explains 
why our solution to the loop equations contains at most terms of order two in 
the fermionic couplings. This is also reminiscent of the bosonic case [25]. At 
the m-th critical point the relevant operators are ctq,<7i, . . .a m -2- The operator 
cr m _i is redundant and the L_i Virasoro condition in the scaling limit can be 
used to compute correlators with one insertion of a m -\. In particular, if in the 
planar limit we evaluate (cr m _i []«)(Zj)) we obtain simply total length L times 
(Yl w (k)) and this loop length counting operator could be used to argue about the 
dependence of planar multiloop correlators on their total length. We believe that 
similar arguments should carry through in our case but using instead the fermionic 
partners of <7 m -i- 

To conclude we should mention that the planar limit of the theory is determined 
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by the following set of equations 



-2k 2 d?F = u 



d t (t + t d t u) 




n 



(5.30) 



n 



n 



n 



Rederiving the one- and higher-point functions from (5.30) is left as an exercise to 
the reader. 



In our approach to the coupling of minimal N — 1 superconformal models to 
2/}-supergravity we have taken as our guiding principle a set of super- Virasoro 
constraints satisfied by the partition function (3.18). This led to our explicit rep- 
resentation of our model in (3.25). From it we derived the superloop equations 
and their planar approximation (3.27), whose solution led to a set of critical points 
labeled by an integer m = 1,2,3, .. . with string susceptibility jst = —1/m and 
with scaling operators a n , with scaling dimensions identical to those expected 
in the NS- and R-sectors of the (2, 4m)-minimal superconformal model. We were 
also able to compute all multiloop correlation functions. There are at least two 
outstanding problems with regard to our model. The first and most important is 
the fact that (3.25) is only an "eigenvalue model" with pairs of eigenvalues (Aj, 9%) 
as basic variables. What kind of generalized matrix model should lead to this 
eigenvalue model is an important open question, and we are quite certain its so- 
lution will not involve supermatrices [26]. The second and easier problem consists 
of proving the fact that the superloop correlators depend only on the total length 
and superlengths by using the properties of redundant fermionic operators. 



6. CONCLUSIONS AND OUTLOOK 
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A physically more interesting question is to find out the properties of super- 
symmetry breaking to all orders in the 1/iV-expansion and non-perturbatively. 
This requires finding the replacement of the Painleve-I equation found in [6] and 
the generalization of Douglas's formulation [7] of the double scaling limit in terms 
presumably of operators realizing a Heisenberg superalgebra. The resolution of 
this problem should also should shed some light on how to generalize our model 
to the coupling of (p,p') -minimal superconformal theories to 2D-supergravity and 
also the explicit form of the c = 1 theory. 

The results presented in this paper lead us to believe that the integrable system 
replacing the KP-hierarchy in our case is not the super-KP-hierarchy constructed 
by Manin and Radul [9] but rather the hierarchy (SKP) found by Rabin [17]. Both 
are defined by starting with a pseudodifferential operator 

L = D + u\ + u 2 D^ + u 3 D~ 2 + ... 

a = .»+{-» (61) 

where we can identify x with the cosmological constant and £ with £1/2- Since 
D 2 = d/dx, D^ 1 = d~ 1 D. Conjugating L if necessary L — > hLh^ 1 in order to 
satisfy Du\ + 2«2 = 0, there exists a unique pseudifferential operator S such that 

S= 1 + siD- 1 + s 2 D- 2 + ... (6.2) 

satisfying 

S^LS = D (6.3) 

The Manin-Radul hierarchy is obtained in terms of a Lax pair constructed in 
terms of L. If the variables t2 n (resp. i2n-i) represent the even (resp. odd) 
flow parameters, and (L n ) + denotes the differential part of L n , the Manin-Radul 



35 



hierarchy is given by 

dL 

dt 2n 
dL 
dt 2n ~ 



1 k=i 



(6.4) 



Note that in the odd flows we have an explicit dependence on the odd parameters. 
In Rabin's case the hierarchy is defined according to the equations 

dS 

(6.5) 



dt 2n 
dS 

dt 2n -i 



-{sd^s- 1 )_s 
-{sd^- 1 s- 1 )_s 



where (L n )- = L n — (L n )+. This hierarchy is integrable but it does not admit 
a representation exclusively in terms of L as a Lax pair. The wave function or 
Baker-Akhiezer function takes the form 



w(z, 9, x, £, t) — z ^exp 



(^~ n + t 2n -i9z- n+1 ) + xz- 1 + £0 



n=l 



(6.6) 



Identifying A with z" 1 and taking 9 to be the same as ours, we see that the 
expression in the exponential is the potential of our model. We believe it should 
be possible to show that our model is equivalent in the continuum limit to this 
hierarchy or some minor modification of it. All these questions are presently under 
investigation. 
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